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We theoretically investigate the applied magnetic field-angle dependence of the fiux-flow resistivity 
pi{ctM) for an uniaxially anisotropic Fermi surface, pt is related to the quasiparticle scattering rate 
7^ inside a vortex core, which reflects the sign change in the superconducting pair potential. We 
find that pf (om) is sensitive to the sign-change in the pair potential and has its maximum when the 
magnetic field is parallel to the gap-node direction. We propose the measurement of the field-angle 
dependent oscillation of pf (qm) as a phase-sensitive field-angle resolved experiment. 
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To elucidate the symmetry of superconducting pair po- 
tential is of great importance for studying Cooper pair- 
ing mechanism in unconventional superconductors (SCs). 
The pair potential is composed of spin and orbital wave 
functions. The orbital wave function is characterized by 
its amplitude and phase (sign of wave function). 

In the last decade, experimental techniques for the 
field- angle resolved specific heat and thermal conductiv- 
ity measurements have developed to identify the Cooper 
pairing symmetry in various superconducting systems 
[ij . These angle- resolved mesurements are powerful tech- 
niques that can detect the anisotropy of the pair potential 
amplitude. However, they cannot probe the sign change 
in the pair potential. That is, these angle-resolved mea- 
surements are not phase-sensitive probes. In addition to 
detecting the anisotropy of the pair potential, it is crucial 
to probe the phase of the Cooper pair in order to discrim- 
inate unconventional SCs including iron-based SCs from 
conventional one. 

Until now, only a few phase-sensitive probes have been 
developed, for example the half-flux quantum observation 
in the tricrystal geometry by scanning SQUID (super- 
conducting quantum interference device) microscope [2| 
and detecting the quasiparticle interference pattern by 
scanning tunneling spectroscopy 3] . In addition to these 
existing experiments, a new phase-sensitive test is highly 
desired. 

In this paper, we propose a new phase-sensitive ex- 
periment, we theoretically study the in-plane field- angle 
dependence of the fiux-fiow resistivity pf (ajvi) for typical 
gap functions and Fermi surface (FS). From our numer- 
ical calculations, we show that the phase-sensitive QP 
scattering inside a vortex core leads to differnt behav- 
ior of pf(aM) between conventional and unconventional 
Cooper pairing. In addition, we show that pf(aM) has 



its maximum when the applied magnetic field is parallel 
to gap-node directions. Our results show that the field- 
angle dependence of the fiux-fiow resistivity can detect 
both the sign change of pair potential and the direction 
of gap nodes. 

The fiux-fiow resistivity pf is dominated by the quasi- 
particle within a vortex core. There are two important 
contributions to pf. One is the QP scattering rate F 
inside a vortex core and the other is the momentum 
dependent quantum level spacing of the vortex bound 
states a;o(fcF) [J|- Here, the QP scattering is due to non- 
magnetic impurities randomly distributed in the system. 

We attribute the fiux-fiow resistivity pf to the energy 
dissipation of the vortex bound states due to the impurity 
scattering inside a vortex core (5| . pf is characterized by 
the two quantities mentioned above [J|, 
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where F^ is the scattering rate in the normal state. We 
have assumed the moderately clean regime 5]. We as- 
sume that the temperature T dependence of pf is pre- 
dominantly comes from 7^ with the QP energy e = k^g^F . 
Here, we have made the rough estimate. Actually the 
QPs distributed with the energy width Ae ^ k^F con- 
tribute to F . The total density of states on a FS is i^o = 
/dS'F/|fF(fcF)| with dSv = |fcF(0fe,^fe)psin6'fcd(/)fe(i6'fc 
being an area element on the FS, the Fermi velocity 
DF(feF) = Vfce(fc)|fc=fcp and the Fermi wave vector fcF = 
|fcF(0fc,6'fe)|(acos0fesin6'fc+6sin(/)fcsin6'fe+ccos6'A:). e(fe) 
is the energy dispersion of electrons. 0fe (0fe) is the az- 
imuthal (polar) angle on the FS. a, h and c denote or- 
thogonal unit vectors spanning crystal coordinates. We 
use the unit system in which ft = 1. 



The momentum dependent inter-level spacing of the 
vortex bound states a;o(A;p) is obtained analytically as 
SQ wo(fcF) = 2|d(fcF)pA^/(|fcF±||vF±(fcF)|) using the 
quasiclassical Green's function method and the Kramer- 
Pesch approximation [8, 9]. We treat the non-magnetic 
impurity scattering by means of the Born approximation 
|10| . The quasiparticle scattering rate for the QPs with 
the energy e inside a vortex core is obtained as fS, Til 
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1 



D{kF,k'F) 



I sine I 



(5) 
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(6) 

Here, (■ ■ • ) = {l/i^o) J cISf ■ ■ ■ /\vFikF)l e(fcF,fc^) = 
6*1, (fep) — dv'ikp) is the scattering angle [see Fig. [T]. T 
has the decay factor exp[— m(soj ^f)] with u(so,fcF) = 

(2|d(fcF)|/|?^F±(fcF)|)/(J'°'ds'^(s')- Fig. n shows the 
quasiparticle trajectories on the plane perpendicular to 
the magnetic field H. The quantities with prime are 
those after scattering. We adopt A{s') = Aoa tanh(s'/^o) 
as the spatial variation of the pair potential, sq and |sq| 
denote the length between the point which is the near- 
est from the vortex core on the QP trajectory and the 
scattering point [3, [ll| . 

In this study, we consider the case that H is applied 
parallel to the a — b plane and rotated. The field angle 
measured from a axis is aM- Here, a, 6 and c are crys- 
tal axes. When calculating the dependence of pf on the 
magnetic field angle aM, we need a_ coordinate system 
fixed to H which is spanned by om, ^m, and cm (vortex 
coordinate ststem) . Here, these axes are orthogonal unit 
vectors with cm set parallel to H (cm || H). vf±, 9y and 
those with prime are defined in the vortex coordinates. 
However, kp and fcp are identified by {(t)k,dk) on a FS 
in the crystal coordinates spanned by a, b and c, which 
characterize the crystal axes. In order to calculate the 
field- angle au dependence of pi, we need to derive the 
relation between vf±, 9y and vf, Ok [1I|- Then, the com- 
ponent of VFikp) projected onto the plane perpendicular 
to H is finally obtained as 

\vF±iA,Ok)\ = \vFi(^k,9k)\n{cj,k,ek), (7) 

^{(t>k,dk) = ycos^flfc +sin^(0fc -aM)sin^6'fc,(8) 
COS e^{(t>k, Ok) = 1 , , „ s| coaOk, (9) 
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Thus, the relation between the vortex coordinate and the 



FIG. 1: The schematic figure of (a) the forward scattering 
and (b) the backward scattering in the vicinity of a vortex, s 
and s' indicate the QP trajectory before and after scattering, 
respectively. 



crystal coordinate is derived. We have now reached to the 
position where we can perform the numerical integration 
on FS and calculate the field-angle dependence of Eq. 

Wc consider the following two types of the simple 
pair potential model. One is line-node s-wave pair: 
d{kF) = I cos2(/)fe| sin^ 0fe. The other is fi^2_y2-wave 
one: d{kF) = cos2(f>kSm^ Ok- Each one has gap nodes 
from the north pole of the FS to the south one in the 
(j)k = (1 + 2n)7r/4 [rad] directions (gap-node directions) 
with n = 0, 1, 2, 3. <pk — ?^7^/2 [rad] directions correspond 
to anti-node directions. In the momentum space, these 
two pair potentials have the same anisotropy in their am- 
plitude |d(fcF)|- However, only the d-wave pair has the 
sign change and the s-wave pair does not. 

Our calculations are performed for an uniaxially 
anisotropic FS with the anisotropy parameter 7 = 
y^mc/mab il2,] . 771c and mab are masses characterizing 
charge transport along the c-axis and within the a ~ b 
plane, respectively. 

We show numerical results for an uniaxially anisotropic 
FS with 7 = 3. In Fig. [2j we show the field- angle depen- 
dence of the fiux-fiow resistivity pi{aM) for the above 
pair potential models. As shown in Fig. [2](a), in the 
case of the line-node s-wave pair, a broad maximum ap- 
pears when H is applied parallel to the gap-node direc- 
tion (ttM = 7i'/4). Note that the field-angle dependence 
of the QP scattering rate /^(ckm) has its minimum when 
H is parallel to the node direction [11|. The oscillation 
amplitude of pf(aM) remains small compared with the d- 
wave case when the temperature T is increased. pf(aM) 
has little dependence on T in the case of the line-node 
s-wave pair. 

On the other hand, in the d-wave case [Fig. HKb)], a 
sharp maximum appears when H is applied to the gap- 
node direction. The oscillation amplitude grows with in- 
creasing T in contrast to the line-node s-wave pair. This 
behavior indicates that the peak of pf (ajvi) has a strong 
temperature dependence in the d-wave case. 

The field-angle dependece of pf is quite contrasting be- 
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FIG. 2: (Color online) Field-angle (qm) dependence of the 
flux-flow resistivity pt in the case of (a) line-node s-wave and 
(b) d-wave pair for a spheroidal FS (7 = 8). Each curve in- 
dicates different temperature. The vertical axis is normalized 
by the minimum value of pf min for each plot. 



tween the line-node s-wave pair and the d-wave one. One 
would question what the reason of this prominent dif- 
ference is. We consider that this difference conies from 
whether there is the sign change in the pair potential or 
not. 

At first, we list the characteristics of /^(fcp)//!!. 
r(feF)//^n is obtained by integrating Eq. (|3]) with respect 
to fcp. Next, we explain the behavior of /^(ckm)- 

In order to clarify why pt^au) behaves contrastingly 
between these pair potential models, we investigate the 
fep dependence of -T first. r{kp)/ F^ indicates which QPs 
are easy to be scattered on the FS. We find the follow- 
ing characteristics of the QP scattering (i) - (iii). (i) 
The QPs in the vicinity of the anti-node direction pre- 
dominantly contribute to F as seen in Figs. [31[a)-(d). 
F{kY)/F^ has higher value around around anti-node di- 
rections (see around dotted lines), (ii) There is the ten- 
dency that the QPs in the direction of H are easy to be 
scattered. As seen in Figs. |3Ua) andEl^b), this property 
of the QP scattering is confirmed by the fact that the 
weight of F shifts a bit toward the field direction. The 
tendency is obvious in Figs. [3][c) andlS^d). The reason 
why QPs have above tendency is because |t;F_L(feF)| of 
QPs in the field direction is small and the contribution 
to F becomes large [see Eq. ([6])]. (in) The other contribu- 
tion to F is expressed by the coherence factor C(fcF, fcp) 
Q, which refiects the sign of the pair potential in Eq. 
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FIG. 3: (Color online) fep dependence of F in the case of (a) 
the line-node s-wave pair when H \\ the gap-node direction, 
(b) d-wave pair when H \ \ the gap-node direction, (c) the line- 
node s-wave pair when H \\ the anti-node direction and (d) 
d-wave pair when H \ \ the anti-node direction on a spheroidal 
FS (7 = 3). The quasiparticle energy e is set to 0.2Zio in 
each plot. The vertical and horizontal axes denote the polar 
angle Qh and the azimuthal angle ^fc respectively. The dotted 
lines indicate the anit-node directions. The field directions 
are indicated by arrows. 
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FIG. 4: (Color online) Field-angle (cvm) dependences of the 
part of F containing the coherence factor and the other part 
of it for d-wave pair with the spheroidal FS. The vertical axis 
is normalized by the minimum value. 



(j4|). This characteristic is discussed in detail in Ref. 
Here we summarize their results in Table HI 



First of all, we consider _r, which is a part of the con- 
tribution to Pi. In the case of the line-node s-wave pair, 
taking into account the characteristics of the QP scatter- 
ing (i) - (iii), we can explain qualitatively the behavior 
of pi(ay\). In the s-wave pair, F(piy\) has its minimum 



TABLE I: The QP scattering types. 



forward scattering backward scattering 



sign-conserved 
sign-reversed 



suppressed 
enhanced 



smaU 
suppressed 



when H is applied to the gap- node direction [ll|. In 
this case, according to the factor (iii), the effect of the 
coherence factor on r{k-p)/rn is small. In addition, in 
the line-node s-wave pair, we confirmed that the coher- 
ence factor has no field-angle dependence. Hence, we can 
neglect the factor (iii) in the s-wave pair. When H is par- 
allel to the gap-node direction, the weight of r{k-p)/ F^ 
shift toward the field direction (i.e., gap- node direction) 
due to the factor (ii). Therefore, the QPs around the gap 
node get easier to be scattered. However, according to 
the factor (i) , the contribution of the QPs in the vicinity 
of the gap nodes to r{k-p)/rn is small. Hence, the large 
contribution to F(fcF)/^n due to the factor (ii) gets small 
due to the factor (i). 

On the other hand, when H is parallel to the anti-node 
direction, the QPs in the direction of H (i.e., anti-node 
direction) , which have small vp±{k-p), can give a large 
contribution to r{kp)/rn due to the factor (ii). In this 
case, differently from the case of JEf || the node direction, 
r remains large due to the factor (i). As a result of 
the above consideration, the minimum of /^(om) appears 
when H is parallel to the gap-node direction. 

In the case of d-wave pair, the backward scatter- 
ing is suppressed [see Table U since the coherence fac- 
tor becomes zero and the scattering point is far from 
the vortex core. However, the forward scattering with 
the sign-change of pair potential is enhanced [Table U 
so that the coherence factor gives the large contribu- 
tion to r{kF)/rn. In Fig. 131 we show the field-angle 
dependences of {{C{kF,k'p)D{kF,k'p))ps')FS, which is 
the part of r{e)/rn containing the coherence factor 
C(A;f, fcp). We also calculate the field-angle dependences 
of ((F(e, A;F,fcp))FS')FS, which does not contain the co- 
herence factor. When H is parallel to the node direction. 



{{C{kF, kp)D{kF, A;p))fs')fs shows the sharp maximum. 
On the other hand, {{F{e, fcp, fcF))FS')FS shows the little 
field-angle dependence. This sharp maximum reproduces 
the behavior of r{aM) in d-wave case \U\ . 

The other contibution to p( is wqCsf)- When 
F(e, fep) =const., the field angle dependence of a;o(feF) 
for a spheroidal FS is qualitatively not different from 
that for a isotropic FS Q. When ujo{kF) =const., 
Pf{T) oc r{e = ksT). In the d-wave case, the dependence 
of wo(A;f) on the field-angle makes the behavior of /^(ajvi) 
stand out. As a result, pf{aM) has a sharp maximum 
when H \ \ the gap-node direction. On the other hand, 
in the s-wave case, the field angle dependence of a;o(feF) 
makes the oscillation amplitude of r{aM) inverted and 
Pf(«M) has its broad maximum when H\\ the gap- node 
direction. 

At last, we comment on the experimental condition for 
measuring the flux-flow resistivity under rotating mag- 
netic field. Our theory is based on the vortex bound 
states formed inside a vortex core. Therefore, extremely 
two dimensional system ,in which a Josephson vortex 
is formed parallel to the layer, is beyond our theo- 
retical framework. In layered organic compounds k- 
(ET)2Cu(NCS)2, in-plane field angular dependence of 
the Josephson-vortex fiow resistance has already mea- 
sured by Yasuzuka et al. |13| . but in three dimensional 
system, the measurement of the fiux-fiow resistivity has 
not performed yet. We calculate Pf (om) also in the case 
of the in-plane anisotropic FS [Tj] . In this case, the be- 
havior of /Of (aivi) is not qualitatively different from that 
in the isotropic FS case. 

In conclusion, we theoretically studied in-plane mag- 
netic field-angle dependence of the fiux-fiow resistivity 
for an uniaxially anisotropic FS. We show that the mea- 
surement of the fiux-fiow resistivity changing the field 
direction within the a — b plane can detect both the 
position of gap nodes and the sign change of pairing po- 
tential. 

The authors thank T. Okada and S. Yasuzuka for help- 
ful discussions. 
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